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Extensional  Vibrations  of  Elastic  Plates' 

p  5 

by  R.  D.  Mindlin  and  M.  A.  Medick^ 


A  system  of  approximate,  two-dimensional  equations  of  extensional 
motion  of  isotropic,  elastic  plates  is  derived.  The  equations  take  into  ^ 
account  the  coupling  between  extensional,  symmetric  thickness-stretch  and 
syianetric  thickness-shear  modes  and  also  include  two  face-shear  modes. 

Hie  spectrum  of  frequencies  for  real,  imaginary  and  complex  wave-numbers 
in  an  infinite  plate  is  explored  in  detail  and  compared  with  the  corre¬ 
sponding  solution  of  the  three-dimensional  equations. 


\ 


^  Htxis  investigation  was  supported  by  the  Office  of  Naval  Research  and 
the  U.  S.  Army  Signal  Engineering  laboratories. 

2  Profe66or  of  Civil  Engineering,  Columbia  University,  New  York,  N.  Y. 

^  Senior  Staff  Scientist,  Avco  Research  and  Advanced  Development  Divi¬ 
sion,  Avco  Manufacturing  Co.,  Lawrence,  Maas.  Formerly  Research  Assistant, 
Department  of  Civil  Engineering  and  Engineering  Mechanics,  Columbia  Uni¬ 
versity,  New  York,  N.  Y. 


2 


Introduction 

In  a  previous  paper  [l]4  a  derivation  was  given  of  approximate 
equations  of  motion  that  take  into  account  the  coupling  between  face  and 
thickness  modes  of  extensional  vibration  of  elastic  plates.  The  thick¬ 
ness  motion  considered  was  that  of  t  .e  lowest  thickness-stretch  mode,  in 
which  the  displacement  is  normal  to  nie  middle  plane  of  the  plate  and 
the  middle  plane  is  the  nodal  plane.  Ihe  frequency  of  this  mode  depends 
on  Poisson's  ratio  (V  ):  the  frequency  Increasing  with  increasing  lA 
When  1/  >  i/J  the  frequency  of  the  lowest  thickness-stretch  mode  is 
hijgier  than  the  frequency,  .Lr dependent  of  IS  ,  of  the  lowest,  symmetric 
thickness-shear  mode.  In  the  latter  the  displacement  is  also  unidirec¬ 
tional  but  it  is  parallel  to  the  middle  plane  of  the  plate  and  there  are 
two  nodal  planes  symmetrically  disposed  with  respect  to  the  middle  plane. 
In  the  range  of  pjisson's  ratios  coranonly  encountered,  both  thickness 
modes  can  couple  with  the  face  modes  and  with  each  other.  This  circum¬ 
stance  has  a  marked  influence  on  phase  and  group  velocities  of  waves  and 
on  the  frequencies  and  shapes  of  high-frequency  vibrational  modes.  In 
the  present  paper  a  derivation  is  given  of  approximate  equations  of 
motion  which  include  the  effects  of  both  thickness-modes. 

The  approximate,  two-dimensional  equations  are  deduced  from  the 
three-dimensional  equations  of  elasticity  by  procedure  based  ox\  the 
series  expansion  methods  of  P>isson  [2]  and  Cauchy  [3l  and  the  integral 
method  of  Kirchhoff  [4].  A  detailed  exposition  of  the  procedure,  using 


4  numbers  in  brackets  refer  to  the  list  of  references  at  the  end  of  the 
paper. 
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a  power  series,  and  its  application  to  approximations  of  orders  zero  and 
one  are  given  elsewhere  [5]  •  In  approximations  of  the  second  and  higher 
orders,  awkward  mathematical  forms  are  encountered  due  to  the  lack  of 
orthogonality  of  the  terms  of  a  power  series.  At  the  suggestion  of  W. 
Prager,  an  expansion  in  a  series  of  Legendre  polynomials  was  studied.  In 
that  case,  although  similar  awkward  forms  appear  as  a  result  of  the  more 
complicated  formula  for  the  derivative,  they  do  not  occur,  for  the  most  J 
part,  until  the  terms  of  the  third  order  are  reached.  Hence,  tire  expan¬ 
sion  in  a  series  of  Legendre  polynomials  is  a  convenient  one  on  which  to 
base  the  present  second  order  approximation.  The  method  of  derivation 
and  the  resulting  equations  (with  the  inertia  terms  omitted)  are  closely 
related  to  E.  Reissner's  theory  [6]  of  three-dimensional  corrections  for 
the  equations  of  generalized  plane  stress. 

Hie  expansion  in  a  series  of  Legendre  polynomials  of  the  thickness- 
coordinate,  followed  by  an  integration  across  the  thickness,  converts  the 
three-dimensional  equations  of  elasticity  into  an  infinite  series  of  two- 
dimensional  equations  which  then  are  truncated  to  produce  the  approximate 
equations.  Die  full  series  expressions  of  displacement,  strain,  stress, 
energies  and  equations  of  motion,  in  conduction  with  an  understanding  of 
Rayleigh's  [7]  exact  solution  of  the  problem  of  vibrations  of  an  infinite 
plate,  are  of  aid  in  deciding  what  to  include  in  various  orders  of  ap¬ 
proximation  and  in  understanding  the  implications  of  what  is  discarded 
and  what  retained.  Hie  series  expressions  and  Rayleigi's  exact  solution 
also  supply  both  the  motivation  and  the  necessary  data  for  making  adjust¬ 
ments,  of  terms  that  are  left  after  the  truncation,  in  addition  to  an 


adjustment  of  the  type  node  by  Bisson  in  establishing  the  zero-order 
aquations.  Bre  adjustment  analogous  to  Poisson's  serves  to  uneouple  the 
nodes  retained  from  the  higher  ones  without  seriously  affeeting  the  be¬ 
havior  of  the  lower  ones.  M  additional  adjustments  are  made  to  improve 
the  match  between  the  frequency  spectra  of  an  infinite  plate  as  obtained 
from  the  approximate  and  emrct  equations.  Brie  is  accomplished  by  the 
introduction  of  coefficients  analogous  to  the  shear  coefficient  in  the 
Timoshenko  beam  equations  [8]  and  the  analogous  equations  for  plates  [91- 
in  the  latter  paper  [9]  it  was  shown  how  the  shear  coefficient  may  be 
chosen  to  effect  a  perfect  match  in  one  or  another  part  of  the  spectrum 
depending  upon  the  frequency  range  and  mode  of  greatest  interest  in  a 
particular  application  of  the  approximate  equations.  In  the  present  case 
four  such  coefficients  are  introduced  and,  due  to  the  complexity  of  the 
spectrum.,  several  possible  combinations  of  matching  points  present  them¬ 
selves.  Bre  choice  is  made  here  to  do  all  the  matching  at  zero  wave  num¬ 
ber.  Bre  range  of  wave  numbers  and  frequencies  over  which  the  match 
remains  good  is  a  measure  of  both  the  usefulness  of  the  approximate  equa¬ 
tions  and  their  range  of  applicability.  In  this  range,  solutions  of  the 
approximate  equations,  in  the  ease  of  finite  plates,  may  be  expected  to 
be  reliable  inasmuch  as  these  solutions  are  composed  essentially  of  com¬ 
binations  of  the  modes  and  overtones  of  the  infinite  plate. 

When  both  the  symmetric  thickness-shear  and  thickness-stretch 
deformations  are  taken  into  account,  important  properties  of  the  fre¬ 
quency  spectrum  contained  in  the  exact  theory  are  reproduced  in  the  re¬ 
sulting  approximate  equations  of  the  second  order,  whereas  they  do  not 
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appear  in  the  previous  approximation  of  the  first  order  [1].  These  prop¬ 
erties  of  the  exact  frequency  spectrum  include  the  imaginary  loop  dis¬ 
covered  by  Aggarval  and  Shaw  [10] ;  the  anomalous  behavior  of  the  second 
and  third  branches  with  variation  of  Poisson's  ratio  [5];  the  frequency 
minimum  of  the  second  branch  [5]  with  its  associated  zero  group- velocity 
at  a  non-zero  wave-number  and  phase  and  group  velocities  of  opposite  sign 
at  smaller  wave-numbers,  as  described  by  Tolstoy  and  Usdin  [11];  and, 
finally,  a  pair  of  complex  branches  which  account  for  edge  vibrations 
observed  in  experiments. 


Expansion  in  Infinite  Series 

We  refer  the  plate  to  rectangular  coordinates  X;  (c  «  1,2,3) 
with  x,  and  x3  in  the  middle  plane  and  the  faces  at  x,  =  *•  b  .  The 

components  of  displacement  Uj  ( j  =  l, 2,3)  are  expressed  as 

where  >j  =Xa/b  ,  the  Pn  (>|)  are  the  Legendre  polynomials 


=  f>  (>|)  *  >|  >  P.  (*l>  *  C3>f*-l)/2,... 

n  ,  .  i  d’MV-O’' 

...  P„(l|)  -  dtj"  ’ 


and  the  Uj<nl  }  it  is  to  be  noted,  are  functions  of  the  coordinates  X,, 

X3  and  the  time,  t  ,  only.  The  are  the  amplitudes  of  polynomial 

distributions  of  displacements  across  the  thickness  of  the  plate.  For 
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convenience,  however,  they  will  be  referred  tv  as  displacements  of  order 
or,  simply,  as  displacements. 

Stress-Equations  of  nation.  The  series  egression  for  the  dis¬ 
placement  is  substituted  in  the  equation 

-°>  (2) 

which  is  obtained  from  the  variational  equation  of  motion  (Reference  [12], 
p.  167).  In  Equation  (2)  the  integration  is  over  the  volume,  V  ,  of  the 
plate;  the  T^j  are  the  components  of  stress;  and  the  summation  conven¬ 
tion  for  repeated  indices  is  employed,  as  are  the  comma  notation  for  dif¬ 
ferentiation  with  respect  to  the  coordinates  X;  and  the  dot  notation 
for  differentiation  with  respect  to  time. 

When  the  integration  with  respect  to  yj  ,  from  - 1  to  +  |  ,  is  per¬ 
formed  in  Equation  (2),  the  result  is 


»Fi"*  -pbC„uJ,"')iuJ<"1dA -0, 


where  A  is  the  area  of  the  plate  and 


Dmn  =  2(n-"r»)  *  1 ,  Cn  =  2/(2n*0  . 
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Hie  t^0'*  and  F- are  defined  as  the  nth-order  coraponents  of 
stress  and  face-traction,  respectively;  while  the  constants  and 
Cn  arise  from  the  operations 


P 

rn-m  j 


O.P.-vfJ"1'"' 

m  =  n  . 


ahe  appearance  of  F.  and  T?.  "  ,  in  Equation  (3),  follows  from 

an  integration,  by  parts,  of  the  terns  in  Equation  (2)  that  contain 
Since  Equation  (3)  must  hold  for  all  A  and  arbitrary  , 

the  quantity  in  parentheses  must  vanish  and  we  arrive  at  the  stress- 
equations  of  motion  of  order  n  : 


bV 


n 

-  z 


i  s-  -  (n) 

=  /obCn  Uj 


In  the  analogous  equations  of  motion  obtained  from  an  expansion  in 
power  series  (Reference  [5],  page  3»04),  an  infinite  series  appears  on 
the  right  hand  side  and  no  series  appears  on  the  left  hand  side.  However, 
the  series  in  Equation  (6)  contributes  more  than  one  term  only  for  to  7  2 . 

Strain.  In  the  three-dimensional  theory,  the  components  of  strain, 
,  are  expressed  in  terms  of  the  components  of  displacement  by 


(7) 


Inserting  the  series  expansion  from  Equation  (l)  and  using  the  formula 

* 

for  the  derivative  from  Equation  (5),  we  find 


e 
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*  !&j  v">p»  *<  v"* 


(8) 


where  <5a.  is  the  Kronecker  symbol  with  i-Z. 

In  order  to  define  components  of  strain  of  order  n  ,  the  summand 
in  Equation  (8)  must  be  expressed  as  the  product  of  Pn  and  a  function 
independent  of  x2  .  Considering  the  double  sura  as  a  triangular  array 
and  interchanging  the  order  of  summation  of  columns  and  rows,  we  find 


(9) 


where  the  ,  defined  as  the  components  of  strain  of  order  n  ,  are 

given  by 


2  4(J"  =  U. U.™  ♦  b"(2„  «n“l (10) 

In  the  analogous  egression  obtained  by  a  power  series  expansion  (Reference 
[5],  P*  3.08),  there  is  no  sura  over  m  .  However,  the  additional  terms  in 
Equation  (..C)  do  not  appear  until  m*  3. 

Stress-Strain  Relations.  The  relations  between  the  T^w'  and  the 
may  be  obtained  by  inserting  the  three-dimensions,],  expressions 


Cn) 


(11) 
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in  the  first  of  Equations  (4). 

T  w  =  Cn 

tj  " 


After  performing  the  integration,  we  find 


(12) 


For  later  use,  it  is  convenient  to  have  Equation  (12)  expressed  in 
the  reduced  indicial  notation,  in  which  double  indices  ranging  from  1  to 
3  are  replaced  by  -ingle  indices,  ranging  from  1  to  6,  as  follows: 


r,  =  f„ , 

t4  -  r„ , 

=  » 

t 

*2  ~  *i1  > 

^3  *  *31  > 

€2  "  (iZ  i 

V 

^  ^31  1 

^3  *  *33  » 

=  Tlt  J 

€3*  €3J> 

v- 

2  . 

Hien  Equation  (12)  becomes 


'  C**q  €q 


in) 


>,v 


6. 


cm  =  cqp  • 


(13) 


Energy  Densities.  Using  the  strain-energy-density,  U  ,  given 
in  the  three-dimensional  theory  by 


2U  -  cijKa  £hJ  =  cp#l  ip  , 
we  define  a  plate-strain-energy-density 

U 


(1«0 


(13) 


and  find,  with  the  aid  of  Equations  (l4),  (9)  and  the  second  of  (5), 


We  also  note  that 


(17) 


(18) 


Similarly,  using  the  kinetic  energy- density,  K  ,  as  given  in  the 
three-dimensional  theory  by 


2K  =  ^UjUj  > 


(19) 


we  define  a  kinetic  energy- density  of  the  plate  by 


—  Jl  -  .  ifi)  •  IfO 

2K  «Zjf(  Kdi,  =  Uj  • 


(20) 


Extensions!  Vibrations  of  Isotropic  Plates 

When  the  plate  is  isotropic,  motions  symmetric  (extensional)  and 
antisymmetric  (flexural),  with  respect  to  the  middle  plane,  may  be  con¬ 
sidered  separately.  In  the  case  of  the  former,  only  those  components  of 
displacement  u j‘n)  are  retained  for  which  )*r\  is  odd.  As  a  re¬ 

sult,  only  those  components  of  stress  T-  1  and  strain  t-*1 
appear  for  which  i  ♦  j  +  o  is  even.  In  the  single  index  notation, 

f  w>  >  €pn)  >  the  terms  which  appear  are  p  +  n  odd  for  p  odd  and 


11 


p*-\n  even  for  p  even. 

°  Hie  stress-equations  of  motion  (6)  then  are 


Hi!!  +  Ik!*  *  f!"  _  2p 

6x,  +  dx3  *  b  "  P  Mr 


bfj*  +  at* 

dX, 


2P 


at1 


ay1  (  ar/"  _  Is!  t  i4<"  _  if  aV™ 

ax,  bxt  b  b  ~  3  fc*s 

(21) 

bf?  rtf'  _  3t£°  2f  &»U,W 

ax,  bx3  b  *  b  5  it1 


at^  b£*  sj*  iT_  gp 
ax,  *  &*3  b  *  b  5 


aC)  at*  st™  tf  .  f  *  .  &  i!naU) 

ax,  *  &x3  to  b  b  ’  7  u‘ 
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Hie  components  of  strain  that  remain,  in  Equations  (10),  ere 


« ;  (u“*  u,'”»  ••■)/b 

=  <su“'/bx,  .  ■£>u,'"/4*i 
«4W  *  »«,”/»«, '•3(u,“tu,"'e-Vb 
tj'1  -  du,"/bx,  »  3(u,“' 

€«  = 

«“  *  J(u“*u,tt>— o/b 

t,"  =  au,7»*. » au.“A«, 

«  ™  a  7(M,W«M,“’e-)/b 

«“  --  au,7K  ■>  7(u,w  *-u,“'*-)/b 


(22) 


ihe  stress-strain  relations  (1.5 )  reduce  to 


r."  *  2[(A<-2/«)e,'“*A(<I'"*<J‘“)] 

t~  -  zr (a •*«,“  *•«,”)] 
t,“  -  »AU,"«0] 

‘V*  1  2/W 
V’  --  4* 

C1  -  2/-  *,7» 

f,al  =  2[(At2/*.)tl“,*A(«,“,*«,")]A 
r,'“  --  2[(A*2/-)f1,u  *A(e,'*.  €,“)]/«■ 
r,“  --  zKA^^le"  ♦  A/(.‘” * 

V”  -- 

t  “  -  2/.  t»/7 

l"'  ‘ZyUt^/7 


where  A  and  yU  Lane's  constants. 
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Hie  strain-energy-density,  in  the  form  given  in  liquation  (17)> 


becomes 


2$  =  r/w</"  +  x ;-61w  ♦  xrm6“ 

*  cv +  cc  ♦cv* 


_  (»j  (t)  r  t*l  OI  ,p  U>  ,  l») 

*Ta  *  l*  *4  +U  *4 


and,  finally,  the  kinetic  energy- density  (20)  is 


K*  p(CirM.w  ^  ca«  ♦  iu^'u,01 


Truncation  of  Series 


We  begin  by  setting 


ur  -  0  ,  u,"  -  (5 ,  n >  2  , 

uaw  «  0,  n>3. 


This  leaves  only  the  components  u,4** ,  Uj'*’ ;  U,1'1,  14,'“,  u,<l) 

and  u.w  ,  as  illustrated  in  Jig.  1.  Hie  terms  u,lm  and  u3‘*' 

are  the  amplitudes  of  uniform  distributions  that  represent  the  thickness 


l-J 


distributions  of  displacements  which  occur  in  low  frequency  extensions! 
and  shear  motions  in  the  plane  of  the  plate;  Ut“'  is  the  amplitude  of 
a  linear  distribution  of  displacement  which  is  an  approximation  to  the 
exact  sinusoidal  distribution  in  the  lowest,  symmetric,  thicxness-stretch 
mode;  U,<M  and  U)"1  are  the  amplitudes  of  quadratic  distributions 
of  displacements  which  are  approximations  to  the  sinusoidal  distributions 
in  the  lowest,  symmetric,  thickness- shear  mode  and  the  face-shear  mode  of 

the  same  order. 

Hie  last  "displacement"  retained  (  Utt%>  )  is  the  amplitude  of  a 
cubic  distribution  which  is  an  approximation  to  the  sinusoidal  distribu¬ 
tion  of  displacement  in  the  second,  symmetric,  thickness-stretch  mode. 
This  is  a  mode  of  higher  order  than  is  to  be  included  in  the  approximate 
equations  of  the  second  order  which  we  seek.  However,  U,'"  produces 
the  second  order  strain  (see  the  eighth  of  Equations  (22))  which, 

in  turn,  appears  in  the  second  order  stress-strain  relations  (see  the 
seventh,  eighth  and  ninth  of  Equations  (23)).  In  order  to  permit  the 
alternating  expansion  and  contraction,  through  the  thickness,  which 
should  accompany  the  stresses  f)<w  and  t}  (by  coupling  through 
FOisson1 s  ratio)  and  at  the  same  time  avoid  coupling  with  the  undesired 
higher  mode,  the  thickness  stress  T/*’  and  the  velocity  u,1  “  are 
set  equal  to  zero.  The  eighth  of  Equations  (23)  is  then  used  to  express 
(.  (t)  in  terms  of  and  (3<l)  : 


(27) 


and  this  result  is  substituted  in  the  expressions  for  and  T 

to  obtain 


r,"  =  2  ve,1*)/*, 

T™  -  ZE'U™  + 


(26) 


where  V  =  \/Z(K+/aa)  is  Poisson's  ratio  and 

t  =>  4/*  (Ay<)/(Afr  Z/a)  »  (29) 

in  which  £  is  Young's  modulus.  In  addition,  the  contributions  of  the 
stresses  7 ft,J  and  Z™  to  the  strain  energy  (see  Equation  (24)) 

are  neglected  in  order  to  destroy  coupling  with  the  unwanted  higher  mode 
when  the  displacements  vary  with  X(  and  . 

At  this  stage,  the  stress-equations  of  motion  end  with  the  fifth 
of  Equations  (21);  the  strain-displacement  relations  end  with  <  w  in 
Equations  (22);  the  stress-strain  relations  end  with  Trm  in  Equations 
(2jj),  with  the  expressions  for  <Tta,>  f2u>  and  f replaced  by  Equa¬ 
tions  (27)  and  (23);  the  strain-energy-density  ends  with  the  third  line  in 
Equation  (24);  and  the  kinetic  energy- density  ends  with  the  fifth  term  in 
Equation  (25). 

Up  to  this  point,  the  process  of  truncating  the  series  expressions 
and  adjusting  the  remaining  terras  is  similar  to  the  one  employed  by 
Poisson  [2]  to  obtain  the  zero-order  equations  of  extensional  motions 
of  isotropic  plates  (Reference  [12],  p.  497);  the  main  difference  being 
that  here  the  process  is  carried  on  at  a  level  two  orders  higher.  In 


1( 


Rjisson's  equations  only  the  zero-order  displacements,  u/*’  and  u,'** 
survive.  These  are  the  amplitudes  of  uniform  distributions,  across  the 
thickness  of  the  plate,  vhich  arc  good  approximations  to  the  nearly  uni¬ 
form  distributions  found  in  the  lowest  extensional  mode  of  the  exact 
theory  (at  long  wave-lengths )  and  are  exact  for  the  lowest  face-shear 
mode  at  all  wave-lengths.  The  additional  terms  Uj10,  u»w, 

which  are  now  included,  are  the  amplitudes  of  first  and  second  degree 
polynomials  and  these  are  not  good  approximations  to  the  distributions 
in  the  thickness- stretch  and  thickness -shear  modes  of  the  exact  theory. 
For  example,  at  infinite  wave-length  along  the  plate,  the  exact  distribu¬ 
tions  are  sinusoidal  across  the  thickness.  It  is  advisable,  therefore, 
to  introduce  additional  adjustments  to  compensate,  as  well  as  possible, 
for  the  omission  of  the  polynomials  of  higher  degrees. 

The  incorrect  distributions  of  displacements  affect  the  frequencies 
mainly  through  the  thickness-strains  and  velocities.  Accordingly,  we  make 
the  substitutions 


.  <•> 

for  , 

», 

tor 

•  in 
wi  1 

ki  *4" 

for  <4  , 

i.  <V" 

for 

u,'” 

for  , 

for 

•  U) 

in  the  strain-energy  and  kinetic  energy  densities,  so  that  the  coeffi¬ 
cients  K»r  (r=  1,2.3,  4)  will  be  available  for  appropriate  adjust¬ 

ments  of  the  equations. 

Finally,  as  a  matter  of  expediency,  we  omit  the  term  utlJ)  from 
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the  strain  .  Ibis  term  complicates  tin  equations  and  may  be  shown 

to  have  little  influence  on  the  long  wave-length  end  of  the  spectrum. 

Second  Order  Approximation 

As  a  result  of  the  truncations  and  adjustments  described  in  the 
preceding  section,  the  variables  and  equations  of  the  second  order  approxi¬ 
mation  are 

Kinetic  Energy- Density : 


+  ov*)] 


m 


St rain- Energy- Density : 


o'"  -  k/cc  *  <,"<") 


t2A(*. <-{,"^*,"0  v-W 


'/“.‘(C'C’CCV3 


(51) 


1 


Stress-Strain  Relations  (  T0ln)  «  dU 


c  = 

2[CA*2/<U/*  r 

A  (ft,  £,- 

*** 

a[»*(A*2/.)6aw' 

♦  Ak,,(*, 

II 

J 

2f  (AtZ^Ur  - 

-  A  (<■/•* « 

II 

% 

*"*  *n 

o  ,  I6> 

V'f 

V*  - 

C  =  <  C/3 

r/*'  -  2 E' («"*♦  V€t'»)/S 
r,"  =  zi'U^^i^Vs 
V*  •-  2/u 

(The  components  of  stress  are  illustrated  in  Fig.  2). 
Strain- Displacement  Relations : 


<•)  , 

•I  “ 

3u, “*/»*, 

tt 

t«>  _ 

* 

u.'Vb 

:  «<  e 

3 

S 

,  <•)  _ 
rJ 

* 

c  ■ 

3u,‘7b  - 

au.'Vsx, 

t  •"  = 

'  3  u,"’/  b 

<“  ■ 

*  »u.'7»«, 

=  au,7>*» 

=  au,7»*. 

-  iu,7ax. 

(The  components  of  strain  are  illustrated  in  Fig.  5). 
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Stress-Equations  of  Motion: 


ar," 

.£* 

--  ip 

dx, 

a  x3 

b 

i&.\ 

at;*’ 

-ip 

ax, 

axs 

to 

aV** 

at* 


a*u3 

at 


<•) 


af, 

ax 


<•» 


at, 

bx. 


<•> 


pW  r  <•/ 

t*  +  -Li¬ 
te  b 


3  at1 


o 


(54) 


ar,'“  »r,“  at)'1  f"  sa**1  jV." 

bx,  ax,  to  b  5"  at* 


bT? 

dx, 


uj 


at, 

ax, 


(j) 


2fi\  h 


tu 


2  P  ^4 

5"  at1 
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Displacement-Equations  of  Motion  (obtained  by  substituting  (33)  in  (32) 
and  the  result  in  (3^)): 


,)^  ♦ 

Xk>, 

W1 

F,m 

-  ■ 

aV“ 

ax, 

b 

ax,  * 

2b 

r 

at* 

Ak,. 

aMt”  , 

Fa'n 

-  o 

axut* 

A  )  -  * 

6x, 

b 

ax, 

Zb 

-  r 

at* 

3A  k,  ec 

3K,*(A  + 

2aOu 

10 

i 

+.  t  i-^L.  -  n  tt3*  J*!ua 

to  2b  r  at* 

.  i-sf/aar+ jur\  tnT 

b  l  ax,  b  i  lb  1  bi 

^\  X, 

>  \  &x,  b  /  ab  r  at* 


(35) 


_ 

b 
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where 


W  dU. 


dx,  *x. 


The  Coefficients  Ku 


The  coefficients  k.>;  are  determined  from  a  comparison  of  prop¬ 
erties  of  the  solutions  of  Equ^tJ-ons  (35),  for  the  case  of  straight- 


crested  waves  in  an  infinite  plate,  with  the  corresponding  properties 
of  the  analogous  solution  of  the  three-dimensional  equations.  Taking 
the  wave-normal  in  the  direction  of  x,  ,  we  find  that  Bjuations  (35) 
separate  into  a  group  of  three  coupled  equations  on  u/**,  U,t0  and 

w  <»>  and  two  independent  equations:  one  on  u,"'  and  one  on  . 

The  coupled  equations  govern  three  "compress ional  modes"  and  the  re¬ 
maining  two  equations  govern  "face-shear"  modes.  To  consider,  first, 
the  group  of  three,  we  set  fj <h)  =  0  and 

U,"»  *  A  sin  $x,  0, 

-  B  cos  fx,  elu,t  ,  (5 

u,w>  -C  s.n^x,  elwt  ,  u3(,,=  0 

in  Equations  (35)  and  obtain  the  secular  equation 


0 
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where 

a„  =  K’i'-n* 

atl  ■-  •  4n,V/ir*  - 

o„  =  fi'/s^A  *  12 «2‘/n*  - 

o„-  2K,(K*-2)z/tT 

«i5  -  -2K./Z/1T 


and 

z=2$b /ir  3  X1*oj/u>,  ,  w,  *  ir{/A/f)n/Zb, 

Kz  *  (A *l/*)//A  =  2 (i-v)/({- 2v)  • 


(38) 


We  also  obtain  the  three  sets  of  amplitude  ratios 


A, 


*i/*i  : 


',2,3 


(39) 


where 


a,  *  t T  (nx-  k  * z^J/2 tc,  (K*-2) 

u  _  _L  /  TTl tC4t.Q>  _  Z'.rfKMjZj* 

‘  ~  2TT\  15  Kt*  /  fn^k1 


and  the  z-  are  the  roots  of  Equation  (37)* 


It  may  be  seen  that  z  ,  which  is  proportional  to  the  wave-number 
f  ,  is  the  ratio  of  the  thicluiess  of  the  plate  to  the  real  half-wave- 
length  along  the  plate;  K  is  the  ratio  of  the  velocity  of  dilatational 
waves  to  the  velocity  of  equivoluninal  waves  in  an  infinite  medium;  and 
fl  is  the  ratio  of  the  frequency  to  the  frequency  of  the  lowest  anti¬ 
symmetric  thickness-shear  mode. 

Hie  relation  between  SI  and  z  ,  in  Equation  (p7),  should  match, 
as  closely  as  possible,  the  corresponding  relation  obtained  from  the  three- 
dimensional  equations.  Hie  match  is  improved,  within  the  framework  of  the 
present  approximation,  by  choosing  appropriate  values  for  the  coefficients 
it,;  ;  but.  since  the  »t;  are  constants,  a  perfect  match  can  be  made  only 
at  one  value  of  Z  for  each  of  them.  Now,  large  enough  *  corresponds 
to  frequencies  high  enough  to  enter  the  range  of  modes  that  have  not  been 
included  in  the  approximate  equations.  In  a  plate  vibrating  at  such  high 
frequencies  under  practical  (i.e.,  not  mixed)  edge  conditions,  the  high 
modes  would,  in  general,  couple  with  the  ones  of  lower  order.  Thus  the 
applicability  of  the  approximate  equations  to  bounded  plates  is  limited 
to  frequencies  below  the  lowest  frequency  of  the  lowest  neglected  mode. 
There  is,  then,  little  advantage  to  be  gained  in  matching  the  approximate 
and  exact  solutions  at  short  wave-lengths  (large  z  )  at  the  expense  of 

9 

a  good  match  at  long  wave-lengths.  In  fact,  we  go  to  the  extreme  and  do 
all  of  the  matching  in  the  neighborhood  of  z  =  0  primarily  because 
of  the  intricate  behavior  of  the  exact  solution  at  long  wave-lengths  and 
also  because  this  choice  results  in  a  reasonably  good  match  out  to  as 
short  wave-lengths  as  the  frequency  limitation  permits. 
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& ten  2  *  0  ,  Equation  (37 )  has  the  three  roots 

rf  •  0 }  J2 K/(*k*/Tr*K,3l  t  6Ck*/ tt1^*  (4o) 

corresponding  to  the  amplitude  ratios 

A  :  B  :  C  =  1:0:0,  0:  l  : 0 ,  0 : 0 1 1 

respectively.  It  may  be  seen  that,  at  infinite  wave-length,  the  three 
"compressional  modes"  uncouple  to  form  an  extensions!.,  a  thickness-stretch 
and  a  thickness-shear  mode.  Ihe  three  limiting  frequencies  are  to  be 
compared  with  the  exact  values 

=  (4i) 

of  the  cut-off  frequencies  of  the  extensions!,  first  thickness  stretch 
and  first  symmetric  thickness  shear  modes  in  Rayleigh's  exact  solution 
[7].  Hence,  we  set 

10, /k3  =  TT/Vil  3  Kt/K4^rr/{is  (42) 

and  these  relations  between  pairs  of  K.-  are  used  in  what  follows. 

We  proceed  to  examine  the  roots  of  Equation  (37)  in  the  neighbor¬ 
hood  of  the  roots  given  in  Equations  (40)  or  (4l)  (which  are  now  the 
same ) . 

First,  for  z  «  I  ?  fi  «  I  ,  we  find 

4z2(k*-l)/k\ 


0O) 


2? 


i.e.,  phase  velocity  f  E/pO-^)]'*  •  This  is  exactly  the ’result  obtained 

fron  the  three-dimensional  theory  for  the  lowest  mode  at  long  wave-lengths. 
Next,  for  Kl,  ,  Equation  (37)  gives 

dfl  _  z  lir'klt1  +  (kz-2)*(k  iil  \  (44a) 

^2  '  kl-4-  i  rile,1  K1  / 


and  in  the  limit,  as  k-*  2  z  —*  0, 


_  +  JUJl i_  .  (45a) 

*  '  "  ^ 


Finally,  for  z  «  i ,  D.  =  2  *  ^  >  £  —*  0  ,  we  find 


dn  _  ir?2  f  4-(kM)(K»-4)  __  J^\  (46a) 

c/z  '  6(k'-4H  5-k*tczl  k.,1  / 


and  in  the  limit ,  as  K  2  and  2  0 , 


dQ  _  _  -n-  #  (47a) 

dz  ' 


Qquations  (44a)  to  (47a)  are  to  he  compared  with  the  exact  values 
(Reference  [5],  p.  2.43)  given  in  the  following  correspondingly  numbered 


(44b) 


equations : 


ii Ls  x  t  ±.  cot^ 

dz  IT  v  4  k  L  ) 
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d  fl  _  ±  2 

dz  ir 


dn 

dz 


z  nr 


TT\Z 


dfl  _  _2_ 
dz  "  +  rr 


(45b) 


(46b ) 


(47b) 


Equating  (44a)  to  (44b)  and  (46a)  to  (46b),  we  find 


ft,1  4«-‘(K -i)(Ka-4)L  k,1  trM2  k  ki 


(43) 


We  also  note,  from  Equations  (43),  that 


1  inn 


192 

rr4 


lino  — 

K-»i  K». 


I  _  3840(TT*-6) 


ir 


(49) 


Equations  (42)  and  (48)  give  the  values  of  the  four  coefficients 
in  terms  of  Poisson's  ratio.  As  a  result  of  these  relations,  all 
of  the  ordinates,  slopes  and  curvatures  of  the  curves  A  vs.  z  , 


characterizing  the  frequency  spectrum  of  the  eompressional  modes  of  the 
approximate  theory,  are  identical  with  those  of  the  exact  solution  at 
Z  0  .  Equations  (48)  and,  in  fact,  all  of  the  equations  of  this 

second-order  approximation  should  be  used  only  for  Poisson's  ratios  less 
than  about  7/l6.  Above  that  value  the  frequency  of  the  thickness-stretch 
mode  is  so  high  in  the  spectrum  that,  in  the  exact  theory,  coupling  with 
the  second  symmetric  thickness-shear  mode  becomes  important  and  that  mode 
is  not  included,  in  the  present  approximation. 

Turning,  now,  to  the  two  face-shear  modes,  we  set  F  —  0 , 

Uj'"  =  A,  cos ,  et“>r)  (50) 

u3a>  =  AjCos$a,  e'wt> 

in  Equations  (35),  and  find  that  u,‘*’  and  u3,a'  are  not,  coupled  and 
have  frequencies  given  by 

nl  *z’,  *  (5D 

respectively.  These  are  to  be  compared  with  the  frequencies 

n'=  z',  4*2.*  (52) 

of  the  face-shear  modes  of  orders  zero  and  two  obtained  irom  the  exact 
equations.  It  may  be  seen  that  the  zero-order  face-shear  mode  is  re¬ 
produced  exactly,  in  the  approximate  equations,  for  all  wave-lengths. 
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The  second  order  face-shear  mode  has  the  correct  ordinate  (A  )  and 
slope  (etfl/dz  )  at  z-0  but  it  does  not  have  the  correct  curvature 
because  of  the  presence  of  the  coefficient  X4  .  By  choosing  K,4  =  / , 
this  discrepancy  could  be  eliminated,  b^t  only  at  the  expense  of  incor¬ 
rect  behavior  of  one  of  the  coupled  modes . 

Frequency  Spectrum 

The  spectrum  of  frequencies  of  the  five  possible  modes  of  vibra¬ 
tion  of  an  infinite  plate,  which  the  approximate  theory  contains,  is 
given  by  the  five  independent  roots  of  Equations  (37)  and  (51)-  We  con¬ 
sider,  first,  the  three  corapressional  modes. 

Equations  (37)  relate  A*  and  x*  but  only  real,  positive 
values  of  A  have  physical  significance.  For  real,  positive  A  ,  the 
roots  xa  may  be  three  real  positive  or  two  real  positive  and  one  real 
negative  or  one  real  positive  and  two  conjugate  complex.  Hence  the  roots 

(55) 

« 

may  be  three  real  or  two  real  and  one  imaginary  or  one  real  and  two  con¬ 
jugate  complex.  In  addition,  the  character  of  the  spectrum  is  different 

according  as  Poisson's  ratio  is  less  than,  equal  to  or  greater  than  l/3. 

VJe  consider  the  case  V1  <1/3  (i.e.,  k  <2  )  first.  Then, 

when  A>  2  ,  the  three  roots  of  Equation  (37)  are  real  (  cf  -  0  ). 

They  are  illustrated  in  Fig.  4a  (for  lA  t/4-  ,  i.e.,  k*  -  3  )  by  the 

three  curves  (full  lines)  marked  <f>t  ,  <f>t  ,  <p 3  in  the  plane  y  =  0 

and  the  region  A  >  2  .  As  A  drops  below  A  •  2.  ,  the  largest 


f 


T 


^9 


root,  ,  approaches  zero  with  asymptotic  behavior  Given  by  Equation 
(4p);  the  intermediate  loot,  ,  approaches  a  minimum,  fl  *  ,  below 
’which  there  are  no  real  roots  other  than  (f> ,  ;  the  smallest  root,  ) 

approaches  a  minimum  in  the  real  plane  0  at  jQ.  *  2  ,  *  “  0 

with  asymptotic  behavior  given  by  Equation  (4oa)  or  (hob)  ’which  are  the 
same  in  •view  of  Equations  ( •  i-O ) .  Continuing  with  the  root  ,  as  _Q 
drops  below  Cl  =  2  ,  the  root  is  imaginary  (  x  -  0  ),  again  with  as- 

yraptotic  behavior  given  by  Equation  (46a)  or  (4bb)  in  the  neighborhood 
of  fl  -  2 .  As  n  approaches  K  ,  from  above,  the  root  ^ 3  for.ns 

a  loop  in  the  imaginary  plane  (  K  -  0  )  and  approaches  _fl 1  k  }  c 0 

with  asymptotic  behavior  given  by  Equation  (44a)  or  (44b).  As  Cl  con¬ 
tinues  to  drop  below  Cl  =  k  ,  the  root  (f)3  becomes  real,  with  behavior 
in  the  neighborhood  of  .Q  »  k  ,  x  -  0  agaii.  given  by  (u4a)  but  now  x 
is  negative.  Upon  further  diminution  of  A  ,  the  root  <p2  approaches 
a  minimum  at  fl-fl*  and  negative  X.  Ihis  portion  of  (p3  (i.e, 

between  fl  =  K‘  ^  and  jQ  =  fl*  )  is  identified  in  Fig.  4a  by  [  <f)j  ], 
where  the  brackets  indicate  that  it  is  the  reflection  in  the  plane  X  *  0 
that  is  shown.  Finally,  when  Cl  <  fl*  3  the  roots  and  <f> 5  are  con¬ 
jugate  complex  (  z  =  ±x  +  iy  ).  One  of  them  is  shown,  in  Fig.  4a,  as  the 
curve  marked  <pz  ;  this  curve  is  also  [  <^3  ],  i.e.,  whe  reflection  of  the 
conjugate  root  <P ^  in  the  plane  x  -  0 . 

As  Poisson's  ratio  approaches  i/p  from  below,  the  frequency  of  the 
thic’ness-stretch  mode  increases,  approaching  the  frequency  of  the  thick¬ 


ness  shear  mode;  i.e.,  at  z.  •  0  the  intercept  Cl  =  k  approaches 
I).  r  2  .  Conjointly,  the  curvature  of  at  O  -  k  }  7.-0 
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approaches  negative  infinity  in  the  plane  i 0  and  positive  infinity 
in  the  plane  *-  0  while  both  slopes  remain  zero:  all  in  accordance 
with  Equation  (44a)  or  (44b).  At  the  same  time,  the  curvature  of  <f>  3 
at  fl=-2  j  z  =  0  approaches  positive  infinity  in  the  plane  -  0 

and  negative  infinity  in  the  plane  0  while  the  slopes  remain  zero: 
all  in  accordance  with  Equation  (46a)  or  (46b).  Meanwhile,  the  imaginary 
loop  shrinks  toward  the  point  H*  2  j  ^  =  0. 

At  "0  -  l/3,  the  thickness-stretch  and  thickness- shear  modes  have 
the  same  frequency  and  the  slopes  of  the  two  branches  of  become  2/  it 

-  in  accordance  with  Equations  (4^a)  and  (47a).  This  situation  is  illus¬ 
trated  in  Fig.  4b.  Here,  again,  the  portions  marked  [  h  ]  are  the  re¬ 
flections,  in  the  plane  X  -  0  of  the  actual  branches. 

When  V  >  1/3  the  thickness-stretch  mode  has  a  frequency  higher 
than  that  of  the  thickness- shear  mode.  The  spectrum  (illustrated  in 
Fig.  4c  for  iX  -  2 /X  )  has  now  undergone  an  important  change  in  that 
the  imaginary  branch  loops  up,  from  XI  ~  instead  of  down. 

Since  only  powers  of  z*  occur  in  Equation  (35 ),  there  is 
another  set  of  physically  significant  roots  given  by  the  reflections  of 
the  curves  of  Fig.  4  in  the  planes  X  =  0  and  ^  -  0 . 

Turning,  now,  to  the  face-shear  modes,  the  first  root  in  Equation 
(51)  yields  the  straight  line  marked  H0  in  Fig.  4.  The  second  face- 
shear  mode  gives  the  roots  marked  Ht  in  the  figures;  these  roots  are 
real  for  Xl  >  2  and  imaginary  for  XI  <  Z  .  As  before,  there  is 
an  additional  set  of  roots  given  by  the  reflections  of  the  curves  Hd 
and  Hj  in  the  planes  X  -  0  and  ^  »  0. 


The  spectra  of  the  corresponding  five  modes,  as  computed  from 
Rayleigh's  solution  of  the  exact  equations,  are  also  shown'*  in  Figs.  4a, 
h  and  c  (as  dashed  lines).  The  importance  of  the  introduction  of  the 
coefficients  *tj  and  their  definitions,  in  terras  of  fbisson's  ratio,  is 
apparent.  Without  these  coefficients,  the  extraordinarily  complicated 
behavior  of  the  branches  of  the  exact  frequency  spectrum  at  long  wave¬ 
lengths  would  not  be  reproduced  in  the  approximate  equations.  In  the 
whole  range  of  frequencies  and  wave  numbers  depicted,  the  approximate 
spectrum  is  reasonably  close  to  the  exact  one ;  the  poorest  representation 
occuring  in  portions  of  the  complex  branches  and  the  spectrum  of  the  sec¬ 
ond  face-shear  mode.  On  the  whole,  fair  results  may  be  expected  from 
solutions  of  the  approximate  equations  in  the  case  of  finite  plates. 

Hie  shapes  of  modes  in  the  various  ranges  of  frequency  can  be 
anticipated  from  the  real,  imaginary  or  complex  character  of  the  roots. 

In  rectangular  coordinates,  for  example,  real  roots  correspond  to  trigo¬ 
nometric  mode-shapes;  imaginary  roots  to  exponential  or  hyperbolic  mode- 
shapes;  and  complex  roots  correspond  to  nodes  whose  shapes  are  given  by 
products  of  trigonometric  and  exponential  or  hyperbolic  functions.  A 
striking  example  of  the  latter  is  to  be  found  in  the  experiments  with 
circular  disks  by  Shaw  [13]. 

The  phase  velocities  (  V  )  and  group  velocities  (  v,  )  can  be 
visualized  readily  from  Figs.  4  inasmuch  as 


^  The  complex  roots  were  kindly  supplied  by  Dr.  Morio  Qnoe. 
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du>  _  djl  ]£7  . 
1  o'* 


(55) 


For  example,  the  phenomenon  of  phase  and  group  velocities  of  opposite 
sign,  noticed  by  Tolstoy  and  Usdin  [11]  in  Rayleigh’s  solution,  is  repre¬ 
sented  by  the  branch  [  3  in  the  real  plane.  Also,  the  minimum  rt 

represents  zero  group  velocity  and  non-zero  phase  velocity.  As  the 
vave- length  approaches  zero  (  x  oo  )  the  frequencies  rise  beyond  the 
range  of  applicability  of  the  equations  and  the  asymptotic  oehavior  of 
the  velocities  of  the  three  compressions!  modes  are 


V!V 


(5*0 


According  to  the  exact  theory,  the  first  of  these  should  be  the  velocity 
of  Rayleigh  surface  waves  and  the  second  and  third  should  be  the  velocity 
of  equivoluminal  waves. 


Additional  Results 

a.  The  equations  of  compatibility  may  be  obtained  by  eliminating 
the  displacements  from  Equations  {5j),  with  the  result: 
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The  first  of  (55)  is  the  usual  compatibility  equation  of  Generalized  plane 
stress  and  the  remaininc  six  equations  correspond  to  the  ordinary  six 
compatibility  equations;  the  main  differences  beinG  that  here  the  compo¬ 
nents  have  a  factor  5  and  the  operator  d/3**  is  replaced  by  l/b. 

b.  It  may  also  be  shown  that  there  are  nine  dislocations  when  the 
strains  and  their  derivatives  are  contir  'ous.  Three  of  them  are  the  two 
translational  and  one  rotational  dislocations  of  Generalized  plane  stress. 
Of  the  remaininG  six,  three  are  translational  (in  the  displacements 
u(“‘,  u/‘,  tO  and  three  are  rotational,  in  the  component  of  rotation 
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c.  A  theorem  of  uniqueness  of  solutions  of  the  approximate 
equations  of  motion  may  be  established  along  the  lines  of  Neumann  s 
theorem  ([11],  p.  176  and  [3]).  Ihis  leads  to  the  following  conditions 
sufficient  for  a  unique  solution  (in  the  absence  of  discontinuities  and 
singularities ) : 

i.  Initial  values  of  U,1*1  ,  ,  U , 1  and 

throughout  the  plate. 

ii.  Initial  values  of  •  u.  and  u, 

throughout  the  plate. 

iii.  One  member  of  each  of  the  five  products  Ftn  u,  ,  F3  w,  * 

F,'"u  ,rt  FWV**,  F  u.“)  at  each  point  of  the  plate,  i.q.,  five  surface 

conditions . 

iv.  At  each  point  on  the  edge  of  the  plate  (normal  m  tan- 
gent  a  )  one  member  of  each  of  the  five  products  l„n  u„  ,  ln,  W,  ,  t„2  u*  , 

i-e->  five  edge  conditions. 

ihe  requirement  that  the  strain-energy-density  \Jl  f  Equation  (31  )> 
be  positive  definite  is  satisfied  by  the  addition  of  the  requirements 
k .,>0;  H,t>  0  to  the  usual  requirements  3 A  +-2/*  >  0  ,  >  0  . 

d.  In  the  case  of  steady  vibrations,  the  displacements  may  be 
expressed  conveniently  in  terms  of  potentials  that  satisfy  Poisson  equa- 

1  ul  t 

tions.  emitting  a  factor  t 


the  results  are 
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/u  V7  *  M0  ♦  Hc  =  0 
yLcV*Ht  +  (n*pujx-  IS b  ) M2  =  0 

VVi  +  ify  '-  «,  i  =  <»2,3 


where 


<*;  *  b  [(A+^m)/-* -pw*-J/*ctA 


(58) 


(59) 


and  the  $t*  are,  again,  the  roots  of  Equation  (57)«  It  may  be  seen 
that  H0  and  Hr  are  the  potentials  of  the  two  face-shear  modes  and 
the  ^>t-  are  the  potentials  of  the  three  compressional  modes. 

e.  Die  tensor,  and  hence  invariant,  characters  of  the  quantities 
that  appear  In  the  second  order  equations  are,  for  the  most  part,  apparent. 


For  example,  ir.  tne  displacement  equations  of  motion  (;p),  if  the  first 
and  second  equations  are  regarded  as  the  rectangular  components  of  a 
vector  and  the  fourth  and  fifth  equations  the  rectangular  components  of 
another  vector,  the  only  differential  operators  that  appear  are  the 
gradient,  divergence,  Laplacian  and  tx/7>tt  :  all  in/ariants.  Ihe 
dependent  variables  are  the  scalar  and  the  two  vectors 


u-.u-k.  +  u-k,  -  u;-'k<  *ur~kr  , 
U‘“  -  u^'K,  ♦u,wk,  •  u;uk,  V  biru,kr  , 
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while  the  gradient  operator  is 


+  *'k  * 


(61) 


where  k,  ,  kA  and  k,  ,  k  „  are  unit  vectors  in  the  rectangular 
directions  1,  5  and  the  orthogonal  curvilinear  directions  cl,  y 
respectively:  all  in  the  plane  of  the  plate. 

The  appropriate  strain  tensors  and  their  expression  in  terms  of 

<* 

vector  displacements  are  not  quite  as  apparent.  We  define  another  vector 
displacement 


(J*  =  t  k,  u, 

<v  <v  ^  i  ► 


tl> 


(62) 


and  a  gradient  operator 


7'=  7  +  k Jt 


(65) 
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where  K  is  a  unit  vector  normal  to  the  plane  of  the  pla^e.  Hien  the 

*4  **  i 

two  tensers 


£w  **(?«“♦  «"?>, 

t'  =  i<v'u'  ♦  u‘V') 


(64) 


constitute  the  strain.  The  seven  equations  of  compatibility,  for  example, 
become 


V  *  * *  7  =  0 

<V  <V  -V  ' 

V'  *  £  '  *v'=0  . 

<v  <v  ^ 


(65) 
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